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Abstract 

The spin-dependent localized Hartree-Fock (SLHF) density-functional approach is extended to the treatment of the inner-shell 
excited-state calculation of open-shell atomic systems. In this approach, the electron spin-orbitals in an electronic configuration 
are obtained by solving Kohn-Sham (KS) equation with SLHF exchange potential and the Slater's diagonal sum rule is used 
to evaluate the multiplet energy of an inner-shell excited state from the single-Slater-determinant energies of the electronic 
configurations involved. This approach together with the correlation potentials and energy functionals proposed by Perdew 
and Wang's (PW) or Lee, Yang, and Parr's (LYP) have been used to calculate the total and excitation energies of inner-shell 
excited states of open-shell atomic systems: Li, B, Ne + , Ne 2+ , Ne 3+ , and Na. The results with the PW and LYP energy 
functionals are in overall good agreement with each other and also with available experimental and other ab initio theoretical 
data. Some new results for highly excited inner-shell states are presented. 

PACS numbers: 31.15.Ew, 32.80.Wr, 32.80.Rm 



I. INTRODUCTION 

Due to computational simplicity and efficiency Q, Q 
density functional theory (DFT) 0, Q has been widely 
applied to many areas in theoretical physics and chem- 
istry as a powerful ab initio approach for the calculation 
of ground-state properties of many-electron systems. The 
basic equation of DFT is Kohn-Sham (KS) equation Q 
and the key part in KS equation is exchange-correlation 
(XC) potential @. Due to incomplete cancellation of 
spurious self-interactions and inherent degeneracy (due 
to the use of spin and angular-momentum independent 
local potentials) of traditional XC potentials obtained 
from uniform electron gas, such as local density approxi- 
mation (LDA) [ll|3| and generalized gradient approxima- 
tion (GGA) Q, Till S||, the conventional DFT using 
LDA or GGA is a ground-state approach. The differ- 
ences of the KS energy eigenvalues of unoccupied and 
occupied orbitals are not rigorously defined as excitation 
energies. However, the KS energy eigenvalues can serve 
as good zeroth-order excited-state energies provided they 
are obtained b y so lving KS equation with a high-quality 
XC potential [l(| • A number of theoretical methods 
have been developed by adopting this point of view E2 • 
In part icular, density work-functional approach (WF) 
[H 1H Q Eg, open-shell localized Hartree-Fock (LHF) 
density- functional approach [l^lL7llL^ll9| . and multiref- 
erence LHF density-functional approach [gjj, etc., 
have been successfully used to calculate excited-state 
properties of atomic and molecular systems. 

Recently, an exchange (X)-only LHF density- 
functional theory has been proposed and successfully ap- 
plied to ground-state calculations of atomic and molecu- 
lar systems E3- I n this X-only DFT, the exchange po- 
tential in the KS equation is a LHF exchange potential 
derived under the assumption that X-only KS determi- 
nant is equal to the Hartree-Fock (HF) determinant. We 



have recently extended this approach to excited states 
of atomic and molecular systems by assuming that the 
X-only KS determinant is also equal to the HF determi- 
nant for excited states [2^]. Based on this postulate we 
have developed a spin-dependent localized Hartree-Fock 
(SLHF) density-functional approach for excited-state cal- 
culation of atomic and molecular systems |22| . In this 
approach, the exchange potential in the KS equation is 
an exact nonvariational SLHF exchange potential con- 
structed for both the ground and excited states. The 
SLHF potential is an analogue of the LHF potential. It 
is self-interaction free and exhibits the correct long-range 
behavior. Further, the SLHF potential requires the use 
of only the occupied orbitals and is dependent of the or- 
bital symmetry of the state. This approach associating 
with Slater's diagonal sum rule [23( has been successfully 
used to calculate multiply excited states of valence elec- 
trons of atomic systems I22J and inner-shell excitation of 
close-shell atomic systems |24J with accurate results. 

In this paper, we extend the SLHF density-functional 
approach to inner-shell excited states of open-shell 
atomic systems. We compute the total and excitation 
energies of inner-shell excited states of open-shell atomic 
systems: Li, B, Ne + , Ne 2+ , Ne 3+ , and Na. In the calcu- 
lation, the correlation potentials and energy functionals 
proposed by Perdew and Wang (PW) [|| and by Lee, 
Yang, and Parr (LYP) || are used to estimate electron 
correlation effect. We will show that the calculated re- 
sults are in overall good agreement with available the- 
oretical and experimental data. We also present some 
new results for the highly excited inner-shell states for 
the first time. 



II. THEORETICAL METHOD 

The SLHF density-functional approach has been dis- 
cussed in Rcf. 221 in detail and is outlined in this section 
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for convenience. 

In spin-dependent density-functional approach, a spin- 
orbital ipi a (r) of the ith electron with spin a (a = a and 
f3 for spin-up and spin-down, respectively) and its orbital 
energy £i a are determined by the KS equation 



H a (Y)ip la (r) = s ia (pia- (r) . 



where, 



(1) 



(2) 



is the KS Hamiltonian and 

Vf (r) = V ext (r) + V H (r) + V xca (r) , (3) 

is the local effective potential. In Eq. V ex t (r) is the 
external potential, Vh (r) is Hartree potential (classical 
Coulomb electrostatic potential between electrons), and 
V XC a (r) is the XC potential. 

For a given atomic system, the external potential 
Vext (r) is known exactly. The Hartree potential Vh (r) 
is given by 



(4) 



where, p (r) = p a (r) + pp (r) is the total electron den- 
sity and p a (r) (for a = a and (3) is the spin-dependent 
electron density defined by 



p° ( r ) = J2 Wi ° ( r )i 2 ■ 



(5) 



Here N a is the number of electrons with spin a and u>i a 
is the occupied number of electrons in the spin-orbital 
<pi* (r). 

The XC potential can be decomposed into the ex- 
change potential V xa (r) and the correlation potential 
V ca (r). In the SLHF density- functional approach, the ex- 
change potential is a SLHF exchange potential V^ sljHF (r) . 
It is given by 



V*r t (r)=V* a (r) + V^(r), 



(6) 



where, 



Vl(r) 



is the Slater potential and 



1 N " f t ct - (r') 



v: 



»,3=1 



(8) 



is a correction to Slater potential. In Eqs. Q and © 

7^ (r) and Qf, are defined by 



7« ( r ) = fia(r)(p j( T(r), 



(9) 



and 



rNL I 



(10) 



where, V^J" is a nonlocal exchange operator of the form 
of HF exchange potential but constructed from KS spin- 
orbitals. 

The SLHF exchange potential determined by Eqs. J|jJ) 
HI has two arbitrary additive constants. The physical 
orbitals can only be obtained by the use of appropriate 
constants in the exchange potential 0. To settle down 
the constants so as to pick up the physical orbitals, it 
is required that the highest-occupied-orbital N a of each 
spin a does not contribute to the correction term V xa (r) . 
In this case, the correction term Vj^(r) decays exponen- 
tially, the SLHF exchange potential behaves asymptoti- 
cally as Slater potential and thus approaches to — 1/r at 
long range 

In atomic systems, an electron spin-orbital is charac- 
terized by three quantum numbers n, /, and a, where 
n and I are the principal quantum number and orbital 
angular momentum quantum number of the electron, re- 
spectively. In the spherical coordinates, the spin-orbital 
ifia- (r) of an electron with quantum numbers n, I, and a 
can be expressed by 



, s Rnla(r) 

<Pia (r) = Y lm ( 



(11) 



where, R n ia{ r ) is the radial spin-orbital, Yi m (0, <f>) is the 
spherical harmonic, m is the azimuthal quantum number, 
and i is a set of quantum numbers apart from spin a of the 
spin-orbital. The radial spin-orbital i? n z<r( r ) is governed 
by radial KS equation, 



2dr 2 



1(1 + 1) 
2r 2 



..off 



('•) 



where i£ ff (r) is the radial effective potential given by 

vf(r) = v ext (r) + v H (r) + „£ HF (r) + v ca (r) . (13) 

In Eq. (H3J), v ext (r), v H (r), v^ HF ( r )> and v ca (r) are 
the radial external potential, radial Hartree potential, 
radial SLHF exchange potential, and radial correlation 
potential, respectively. 

For an atomic system with a nuclear charge Z, the ex- 
ternal potential is the Coulomb potential between elec- 
tron and nucleus 



Vext (r) 



(14) 



In central-field approach, the radial Hartree potential 
is calculated from 



v H (r) = 47r / —p(r')r' 2 dr', 



(15) 



where, r> is the larger of r and r', p(r) = p a (r) + pp(r) is 
the spherically averaged total electron density, and p a (r) 
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(a = a or 0) is the spherically averaged spin-dependent 
electron density given by 



1 

47T 



p a {v)dVt 



1 " 

— y 

nl 



W n la 



Rm 



(16) 



Here the symbol v a stands for a set of quantum num- 
bers for summation and the sum is performed over all 
the occupied spin-orbitals with spin a. This expression 
is accurate for spherically symmetric (close-shell) states, 
but it is only an approximation for non-spherically sym- 
metric (open-shell) states. It may induce an error when 
it is used to evaluate the energy of a non-spherically sym- 
metric state. However, the error is negligible compared 
to the order of calculated multiplet splitting 11]. 
The radial SLHF exchange potential is given by 



,,SLHF 



(r) = *4 ( r ) + v c a (r ) ) 



where, 



vL(r) 



ATtp a (r) 



EE 

nlm n'l'm' 



(17) 



(18) 



is the radial Slater potential and 



(r) 



E 



nlm n I'm 



ul m ,n' V in 



,(r). 



(19) 



is a correction to the radial Slater potential. The matrix 
elements s£ lm _», m ,(r) and c a nlm n , Vm ,(r) in Eq. (SJ are 
given in Ref. |22j. 

To calculate electron spin-orbital, the Legendre gen- 
eralized pseudospectral (LGPS) method [25| is used to 
discretize the radial KS equation (|12Jl . This method asso- 
ciated with an appropriate mapping technique can over- 
come difficulties due to singularity at r = and long-tail 
at large r of Coulomb interaction and thus provides a 
very effective and efficient numerical algorithm for high- 
precision solution of KS equation. Using the electron 
spin-orbitals of an electronic configuration, a single Slater 
determinant for a specific electronic state is constructed 
and its total energy calculated. The total energy is a 
sum of non-interacting kinetic-energy Ek, external-field 
energy E ext , Hartree energy Eh, exchange energy E x , 
and correlation energy E c . The values of Ek, E ext , Eh, 
and E x are evaluated by 



Rnia (r) dr, 



a— a nl 

1(1 + 1] 



~2dr 2 



2r 2 



Eact = 47T 



tt (r)p(r) r 2 dr, 



(20) 



(21) 
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(22) 



(23) 



where, n represents a collection of all the quantum num- 
bers involved, the matrix elements n?L v^i, F„ k 



and G ilv,n'l'a 



Vlnifl'Tri'i ± nla,n'l'< 

are given in Ref. j2^ ■ 
For a multiplet state that can be described completely 
by a single Slater determinant, the energy is calculated 
directly from the single Slater determinant. For a multi- 
plet state that cannot be represented by a single determi- 
nant, the energy can be calculated by means of Slater's 
diagonal sum rule j2^| • According to this rule, a sum over 
single-Slater-determinant energy E{Y>i) of determinant 
Dj from an electron configuration equals to a weighted 
sum over multiplet energy E(Mj) of multiplet state Mj 
involved in the same electron configuration, namely, 



(24) 



where, the weight dj is the times that the multiplet state 
Mj appears in all the single Slater determinants. Similar 
procedures have been employed in recent excited-state 
calculations Hill, El. 



III. RESULTS AND DISCUSSION 

The procedure described in the preceding section is ex- 
tended to calculate the total energies (E) and excitation 
energies (Ai?) of inner-shell excited states of open-shell 
atomic systems: Li, B, Ne + , Ne 2+ , Ne 3+ , and Na. In the 
calculations, the correlation effect, which is character- 
ized by the correlation potential v ca (r) and correlation 
energy E c , is taken into account through the correlation 
potentials and energy functionals of Perdew and Wang 
(PW) 9] and of Lee, Yang, and Parr (LYP) 0, respec- 
tively. The results obtained with these two correlation 
energy functionals are listed in columns PW and LYP in 
the following tables, respectively. 



A. Inner-shell excitation of Li 

Inner-shell excitation of Li has been the subject of ex- 
tensive experimental and theoretical studies. Both exper- 
imental [p, jp, jp, JM Jpjpjpj3lj3|| and theoretical 

data IlilalllliMMliM 

are abundant. In TABLE [Q to TABLE |Hl] we present 
the total and excitation energies of our calculations with 
PW and LYP correlation potentials and energy function- 
als for inner-shell excited states ls2sns 2,4 S (n = 2 ~ 8), 
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Is2snp 2 > 4 P (n = 2 ~ 8), and ls2p 2 2 S, 2 < 4 P, and 
2 D, respectively. For comparison we also list in these 
tables representative experimental data (2^, EI EH as 
well as theoretical results obtained from calculations us- 
ing density work- functional formalism (WF) |36l , saddle- 
point complex-rotation (SPCR) [37l Isslls^ l4pj . extensive 
configuration-interaction calculation (ECI) [43, 144 [iM 
146] , and Rayleigh-Ritz variational method f RRV)|49| . 

For the total energies of inner-shell excited states 
ls2sns 2 ' 4 S listed in TABLE HJ the maximum relative 
discrepancies of our PW and LYP results are 0.10% and 
0.32% to the WF results, 0.21% and 0.43% to the SPCR 
results, and 0.19% and 0.62% to the ECI results, respec- 
tively. This indicates that both the PW and LYP re- 
sults are in good agreement with WF, SPCR, and ECI 
results, and the PW results are a little bit better than 
the LYP results. For the excitation energy of excited 
state ls2s 2 2 S, the maximum relative deviations of our 
PW and LYP results and the WF results are 0.91%, 
1.65%, and 1.25%, respectively, to the experimental re- 
sults. This demonstrates that the PW result is better 
than both the LYP and WF results. For excited states 
ls2sns 2 S (n > 3), other theoretical results of excitation 
energy are not available. For these states the maximum 
relative discrepancies of our PW and LYP results of ex- 
citation energies to the experimental results are 1.45% 
and 2.75%, respectively. For the excited state Is2s3s 4 S, 
the maximum relative deviations of our PW and LYP 
results of excitation energies to the experimental result 
are 0.75% and 1.60%, respectively. While the maximum 
relative discrepancies of the ECI and RRV results to the 
experimental result are 0.59% and 0.50%, respectively. 
Hence our PW result is in good agreement with both 
the experimental result and other theoretical results for 
this state. For excited states ls2sns A S (n > 3) there 
is no experimental data available. For these states, the 
maximum relative deviations of our PW and LYP results 
to the RRV results are 0.68% and 1.39%, respectively. 
Thus our PW results are close to the RRV results. For 
the all the results above the PW results are better than 
the LYP results. As will be shown, due to overestima- 
tion of LYP energy functional to correlation energies of 
atomic systems with smaller Z || , the LYP results 
are generally worse than the PW results for atomic sys- 
tems with Z < 8. 

For the total energies of excited states ls2snp 2 ' 4 P 
given in TABLE ITT1 the maximum relative discrepancies 
of our PW and LYP results are 0.56% and 0.89% to 
the WF results, 0.64% and 1.03% to the SPCR results, 
and 0.25% and 0.82% to the ECI results, respectively. 
This illustrates that both the PW and LYP results are 
in agreement with the complicated ab initio methods. 
For the excitation energy of excited state \s2s2p 4 P, the 
maximum relative deviations of our PW and LYP re- 
sults to the experimental results are 0.38% and 1.66%, 
while the maximum relative discrepancies of the WF, 
ECI, and RRV results to the experimental results are 
1.423%, 0.03%, and 0.07%, respectively. For excitation 



energy of the excited state \s2s2p 2 P, the maximum rel- 
ative discrepancies of the PW, LYP, and WF results to 
the experimental result are 0.05%, 1.70%, and 1.85%, re- 
spectively. For excitation energy of the state Is2s3p 2 P, 
the maximum relative deviations of the PW, LYP, and 
WF results to the experimental result are 1.35%, 2.73%, 
and 1.11%, respectively. Thus the PW results are better 
than the LYP results and a little bit better than the WF 
results. 

For the total energies of inner-shell doubly excited 
states \s2p 2 2 S, 2 - 4 P, and 2 D given in TABLE D]]] the 
maximum relative deviations of our PW and LYP re- 
sults are 0.63% and 0.75% to the WF results, 0.47% and 
0.78% to both the SPCR and ECI results, respectively. 
This demonstrates that the PW energy functional has al- 
most the same precision as the LYP energy functional in 
calculation of total energies of these inner-shell excited 
states. For excitation energies, the relative discrepancies 
of our PW and LYP results to the experimental result 
are less than 1.10% and 1.50%, respectively, while the 
maximum relative discrepancies of the WF results to the 
experimental results are 1.09%. Hence the PW results 
are very close to the WF results and a little bit better 
than the LYP results. 



B. Inner-shell excitation of B 

The second open-shell atom for which the total ener- 
gies of inner-shell excited states are computed is B. In 
TABLE IIVI we present Auger transition energies from 
the inner- shell excited states ls2s 2 2p 2 2 S, 2 < 4 P, and 2 D 
of B to the singly excited states ls 2 2s2p *' 3 P and the 
doubly excited states ls 2 2p 2 1 S, 3 P, and l D of B+. For 
comparison we als o g ive in this table the available ex- 
perimental results |34| as well as theoretical results from 
WF [111, perturbation Z-expansion theory (PZE) jMEl, 
and perturbation Z-expansion theory with relativistic ef- 
fects (RPZE) HJ. Except for the transitions ls2s 2 2p 2 
2 D -> ls 2 2p 2 l S and ls2s 2 2p 2 2 S ls 2 2p 2 1 S, the 
maximum deviations of our PW and LYP results to the 
experimental results are 0.634% and 0.241%, while the 
maximum discrepancies of the WF, PZE, and RPZE re- 
sults to the experimental results are 0.783%, 1.211%, 
and 0.4091%, respectively. This indicates that our re- 
sults are much better than the PZE results, a little bit 
better than the WF results, and much closer to the 
RPZE results. The maximum deviations of the PW, 
LYP, PZE, and RPZE results to the experimental results 
are 1.821%, 1.033%, 3.257%, and 0.07% for the transition 
ls2s 2 2p 2 2 D \s 2 2p 2 X S and 0.437%, 1.232%, 1.392%, 
and 0.04% for the transition ls2s 2 2p 2 2 S -»• \s 2 2p 2 X S, 
respectively. This indicates again that our results are 
better than the PZE results. In addition, according to 
our result the spectrum line with 170.7 eV previously 
signed to the transition ls2s 2 2p 2 2 D — > \s 2 2p 2 X S may 
likely belong to the transition \s2s 2 2p 2 2 D — > ls 2 2p 2 3 P 
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TABLE I: Total energies (E) and excitation energies (AE) of inner-shell excited states ls2sns 2 ' 4 S (n — 2 ~ 8) of Li. The 
ground state energies obtained from calculations with PW and LYP correlation potentials and energy functionals are —7.4837 
(a.u.) and —7.4872 (a.u.), respectively. Here 1 a.u. =27. 2116 eV is used. 
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TABLE II: Total energies (E) and excitation energies (AE) of inner-shell excited states ls2snp 2,4 P (n = 2 ~ 8) of Li. 
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C. Inner-shell excitation of positive ions Ne + , 
Ne 2+ , and Ne+ 

To explore the feasibility of the approach to inner-shell 
excitation of atomic ions, we also apply the procedure 
to inner-shell excited-state calculations of positive ions 
Ne+, Ne 2+ , and Ne 3+ . In TABLE |V| we present the 
excitation energies of optically allowed transitions in- 
volved in inner-shell excited states of these ions along 
with the theoretical results of multiconfiguration Dirac- 
Fock (MCDF) method [12 . Note that the MCDF results 
here are the weighted-averaged values of those given in 
|52| . The maximum relative deviations of our PW and 
LYP results to the MCDF results are 0.171% and 0.071% 
for Ne+, 0.162% and 0.059% for Ne 2+ , and 0.159% and 
0.081% for Ne 3+ , respectively. Thus the agreement of 
our results with the MCDF results are quite satisfac- 
tory, demonstrating that the SLHF together with PW 
and LYP correlation potentials is accurate for the calcu- 
lation of inner-shell excited states of Ne positive ions. In 
addition, the LYP results are more accurate than the PW 



results for these atomic systems with Z = 10 p2l |2^|. 

D. Inner-shell excitation of Na 

Finally we compute the total energies and excitation 
energies for inner-shell excited states ls2s 2 2p 6 3s( 1 S)np 
2 P ( n = 3 ~ 8) and ls2s 2 2p 6 3s( 3 S)np 2 ' 4 P (n = 3 ~ 8) 
of Na. The results are shown in TABLE together 
with the theoretical results from relativistic configura- 
tion interaction method (RCI) [5j| and multiconfigura- 
tion Hartree-Fock method (MCHF) ,54] as well as the 
experimental results (Exp.) j55j for comparison. 

For the excitation energies of excited states 
ls2s 2 2p 6 3s( 1 S)np 2 P, the relative deviations of 
our PW and LYP results to the experimental results 
are not more than 0.18% and 0.07%, respectively. This 
indicates that both the PW and LYP results are in good 
agreement with the experimental results and the LYP 
results are better than the PW results. In addition, the 
maximum relative discrepancies of our RCI and MCHF 
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TABLE III: Total energies (E) and excitation energies (AE) of inner-shell doubly excited states ls2p 2 2 S, 2 ' 4 P, and 2 D of Li. 
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"PW results, *LYP results, c 36], d 38, 40], " 43, 44, 45], f 49], 9 29], * 31], and * 34]. 



TABLE IV: Auger transition energies (AE) from inner-shell excited states ls2s 2 2p 2 2 S, 2 ' 4 P, and 2 D of B to singly excited 
states ls 2 2s2p 1,3 P and doubly excited states ls 2 2p 2 1 S, 3 P, and 1 D of B + . 



Initial states 


Final states 






AE 








(eV) 




Present 




Other 






theory 


jlIj xpci unc hi 


PW" 


LYP° 


WF C 


PZE a 


RPZE° 


RBB J 


ls2s' A 2p A 


4p 


\s 2 2p 2 3 


P 


168. 


3500 


169.1008 




167 


.4101 


169.10 


169.2 






ls 2 2p 2 1 


D 


167. 


.3105 


168.5612 




166 


.7952 










ls 2 2p 2 1 


S 


165. 


5064 


167.2558 




163 


.2454 








2 D 


\s 2 2p 2 3 


P 


170 


.4344 


170.7805 




169 


.3045 










\s 2 2p 2 1 


D 


169 


3949 


170.2409 




IBS 


.6895 










\s 2 2p 2 1 


S 


167. 


5908 


168.9356 




165 


.1397 


170.58 


170.7 




2p 


\s 2 2p 2 3 


P 


170 


.4181 


170.2883 


171.5858 


169 


.0290 


170.73 


170.7 






\s 2 2p 2 1 


D 


169. 


3786 


169.7487 




168 


.4140 










\s 2 2p 2 1 


S 


167 


5745 


168.4433 




164 


.8643 








2 S 


\s 2 2p 2 3 


P 


172 


.7828 


173.1289 




171 


.0097 










\s 2 2p 2 1 


D 


171. 


.7433 


172.5893 




170 


.3947 










\s 2 2p 2 1 


S 


169 


9392 


171.2840 




166 


.8450 


169.13 


169.2 


ls2s 2 2p 2 


4 P 


ls 2 2s2p 


3p 


176. 


1516 


176.8376 




175 


.7192 










ls 2 2s2p 


X P 


171. 


4249 


172.9961 




170 


.8674 








2 D 


ls 2 2s2p 


3p 


178. 


2360 


178.5173 




177 


.6135 










ls 2 2s2p 


X P 


173. 


5093 


174.6759 


173.2321 


172 


.7617 


174.17 


174.6 




2p 


ls 2 2s2p 


3p 


178. 


.2197 


178.0251 




177 


.3381 










ls 2 2s2p 


l P 


173 


.4930 


174.1836 


173.4062 


172 


.4862 


173.90 


174.6 




2 S 


ls 2 2s2p 


3p 


180 


.5843 


180.8657 




179 


.3188 










ls 2 2s2p 


1 P 


175. 


8577 


177.0242 




174 


.4669 
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TABLE V: Excitation energies (Ai5) for optically allowed 
transitions of inner-shell excited states of Ne fc+ (k = 1 ~ 3). 
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"PW results, b LYP results, and c obtained with data in 1521 . 



results to the experimental results are 0.06% and 0.04%, 
respectively. Thus our LYP results are very close to the 
RCI results and MCHF results. Furthermore, according 
to our results the identification of the photoionization 
spectra ls2s 2 2p 6 3sC_S)3p 2 P and ls2s 2 2p 6 3s( 3 S)3p 2 P 
in the experiment £3 should be exchanged. This result 
agrees well with the theoretical results of [Hi) and ■ 
For the excitation energy of excited state 



ls2s 2 2p 6 3s( 3 S)3p 2 P, the maximum relative devia- 
tions of our PW and LYP results to the experimental 
results are 0.22% and 0.16%, respectively, which are a 
little bit larger than those of RCI and MCHF results to 
the experimental results. For the excitation energies of 
excited states ls2s 2 2p 6 3s( 3 S)np 2 P with n = 4 and 5, 
the maximum relative deviations of our PW and LYP 
results are 0.12% and 0.04% to the RCI results, and 
0.18% and 0.10% to the MCHF results, respectively. 
This demonstrates that our LYP results are very close 
to the RCI results for these states. 



For the excitation energies of excited states 
ls2s 2 2p 6 3s( 3 S)np 4 P, the relative deviations of our PW 
and LYP results to the RCI results are less than 0.13% 
and 0.03%, respectively. This illustrates again that 
our LYP results are in very agreement with the RCI 
results for these states. For all the excited states the 
LYP results are better than the PW results. Thus the 
LYP energy functional is more accurate than the PW 
energy functional in the excitation energy calculation of 
inner-shell excited states of atomic systems with large 
Z. 
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TABLE VI: Total energies (E) and excitation energies (A_E) of inner-shell excited states ls2s 2 2p 6 3s( 1 S)np 2 P and 
ls2s 2 2p 6 3s( 3 S)np 2A P (n = 3 ~ 8) of Na. The ground state energies obtained from calculations with PW and LYP cor- 
relation potentials and energy functionals are —162.2265 (a.u.) and —162.2687 (a.u.), respectively. 
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IV. CONCLUSIONS 

In summary, the procedure developed for excited-state 
calculation based on the SLHF density functional ap- 
proach and Slater's diagonal sum rule has been extended 
to the treatment of inner-shell excited states of atomic 
systems. In this procedure, the electron spin-orbitals in 
an electronic configuration are obtained first by solving 
the KS equation with the exact SLHF exchange potential. 
Then a single-Slater-determinant energy of the electronic 
configuration is calculated by using these electron spin- 
orbitals. Finally, a multiplct energy of an excited state is 
evaluated from the single-Slater-determinant energies of 
the electronic configurations involved in terms of Slater's 
diagonal sum rule. The key part of this procedure is the 
SLHF exchange potential. We have applied this proce- 
dure to the calculations of total energies and excitation 
energies of inner-shell excited states of open-shell atomic 
systems: Li, B, Ne + , Ne 2+ , Ne 3+ , and Na. The corre- 



lation effect is taken care of by incorporating the PW 
and LYP correlation potentials and energy functionals 
into calculation. The results from the calculations with 
LYP and PW energy functionals and energy functionals 
are in overall good agreement with each other and also 
with the available more sophisticated ab initio theoret- 
ical results and experimental data. This demonstrates 
that the SLHF density-functional approach can provide 
a simple and computationally efficient approach for the 
accurate calculation of inner-shell excited states of open- 
shell atomic systems within DFT. 
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